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Abstract 

Using the recently amended sea-boson method, we compute the mo- 
mentum distribution of the one-band Hubbard model in one and two 
spatial dimensions. We compute the asymptotic features of the momen- 
tum distribution explicitly away from half filling for weak coupling in one 
and two dimensions. While the results are not exact by any means, they 
provide the exact asymptotics, namely they are able to reproduce the ex- 
ponents obtained by Shulz in one dimension obtained using Bethe ansatz. 
The corresponding results in more than one dimension are therefore as 
believeable. 



1 Introduction 

In the present article, we use the recently perfected sea-boson method [1] [2] [3] 
to compute the momentum distribution of the one-band Hubbard model. The 
Hubbard model in two dimensions is thought to be important in understanding 
high-T c superconductivity. Thus a proper solution of this using the sea-boson 
method can shed new light on the question of breakdown of Fermi liquid be- 
haviour in two dimensions from which one may decide whether or not the Hub- 
bard model is adequate in describing high-T c materials. Furthermore, lattice 
Fermi gases are easy to solve using the sea-boson approach since the hamilto- 
nian remains separable even after invoking the repulsion attraction duality [2] [3]. 
The anomalous exponent is computed for the Id case and a comparison is made 
with the exact Bethe ansatz results of Lieb and Wu[4] as elaborated by Shulz [5]. 
Unfortnately as in the usual bosonization in Id [7] we are unable to probe the 
strongly coupled ( large U ) Hubbard model since the velocities of the spinons 
becomes imaginary and the formalism breaks down. However we take the plau- 
sible point of view, supported by earlier works [3] , that says that breakdown 
of Fermi liquid theory (in any dimension) comes from the infrared divergence 
of certain integrals in the sea-boson formalism and not from the strength of 
the coupling. Thus it is sufficient to investigate the weakly coupled regime to 
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ascertain whether or not Fermi liquid theory has broken down. This point of 
view is further strengthened by the remarks of Shulz[5]. 

The Hubbard model is one of the most extensively studied models in Con- 
densed Matter Physics. The collection of reprints in the volume by Korepin and 
Esslcr[8] is particularly useful. Also the monumental works by Weng ct.al.[ll] 
which is a path integral approach to the Id Hubbard model cannot be ignored. 
In more than one dimension, there are many attempts notable among them are 
by Metzner[9] and Fukuyama [10]. 



2 Hubbard Model in One Dimension 

Here we write down the one-band Hubbard model in one dimension in momen- 
tum space. To do this we must be careful since there are issues of backward 
scattering and umklapp process that play an important role the latter especially 
at half-filling[6]. Consider the model written in the usual manner. We assume 
that we are in the thermodynamic limit so that we may assume that the number 
of sites N a 3> 1 is even. Also we set the hopping matrix element t = 1 and we 
work with units such that the lattice spacing a = 1. This means L — N a . 

m=N a /2-l N a /2-l 
H = - E E ( C lna c m+la + c] n+1(J C ma ) + U E i 1 ) 
m=-N a /2 cr m=-N a /2 



Here Tima <^ rna c raa . 



c ma = -±=J2e ikm c ka (2) 
VJVa k 

JVa/2-l 

C fecr = prj- ^ ' e C mcr (3) 



v ^ a m=-N a /2 

Since we want periodic boundary conditions we must have, 

Cm+N a a — C mcr (4) 

This means that k ■ N a = 2ir v where v = 0,±1,±2, We can also see 

that Ck+2TT,a = This means that we may now assume that in the ther- 

modynamic limit, the k's form a continuum and take on values in the interval 
k G (— 7T, 7r) = /. We may do so without any loss of generality. This perdiodic- 
ity leads to some novel processes such as umklapp processes. To see this more 
clearly we write down the repulsion term in momentum space. The hopping 
term is straightforward. To see this we write, 

H kin = E E e * c L c ^/( fc - fc ') + H - c - ( 5 ) 

° k,k' 
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where, 



/ ( »-V)--» ^y> -. 



if fc - k ^ 0, ±2tt, ±4tt, ... and 

f(k-k) = N a (7) 

if fc — k' = 0, ±27r, ±47r, .... Since k, k' e (— w, it) there is no chance for k — k' = 
±2ir. Thus k — k = is the only possibility. However in the repulsion term 
we find that such umklapp processes are indeed possible. The hopping term is 
therefore given by 

H k in = ^2 e k c\ a c k a (8) 
fee/, ci- 
and efe = —2 cos(k). The interaction term is given by, 

m=N a /2-l 

H v = U ^ n m^n. mi (9) 

m=-N a /2 

In momentum space it is, 

H " = ^2 EE c i'T Cfe t c J/i c ^ /(fc - k ' +p-p) ( 10 ) 

k,k p,p 

N a /2-l 

f(k-k +p-p)= e lm(k - k+p - p) (11) 

m= — N a /2 

We can see that unless p = 0, ±2ir, ±4ir, .... we must have f(p) = 0. Thus we 
may write, 

f(k-k' +p-p) = N a (5 k _ k , +p _ p , a + 5 k _ kl+p _ p/27I + 6 k _ k , +p _ p ,_ 27r ) (12) 

Since k,k ,p,p E (— ir, it), there is no chance for k—k +p—p = ±4ir or anything 
higher ( since — 4tt < k — k + p — p <47r). However, k — k + p — p = ±2w is 
perfectly possible, and this is precisely the umklapp process. Therefore we may 
write, 

Hu = C i'| Cfe T C p'^ c Pl(^fc-fc'+p-p',0 + ^k-k'+p-p' ,2tt + ^k-k'+p-p ,-2tt) 

k,k' p,p' 

(13) 

In other words, 
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+ W a E C l+ q ^ c l- q +Gl C Pl ( 14 ) 

k,p,q 

The last two terms correspond to umklapp processes. Here G = 2n is the 
reciprocal lattice vector. Define the sea-displacement operator, 



pa 



n F (p)(l-n F ( P ')) cl (15) 



= c]p a c piJ and Uf{p) = 0(cf — £p)- Eq.( 15) may be formally inverted 
and a formula for the number conserving product of two Fermi fields may be 
written down. In the RPA-sense we may regard the object A** as being 
small in the sense that the matrix respresentation of this operator in a suitably 
restricted Hilbert space is sparse. Thus we can be content at including only the 
leading terms. The formula for the off-diagonal product c| CT c p ' a > for p ^ p may 
be written down as follows. 

4*<Vo"»4f'+ A %> ( 16 ) 

If p = p we have instead 1 , 

c\ aCpa , = 5^,n F {p) + J2 4 P -V^V - E Jp^A&v* (17) 
qoi q<y\ 

The interaction term with q = may be omitted since this leads to a c-number. 
The terms with k = k ± q are also omitted since they result in a hamiltonian 
that is quartic in the sea-bosons (see later for some problems related to this). 
Since these objects are small in the restricted Hilbert space we may ignore these 
terms. The full hamiltonian may be written as the sum of two terms. The first 
corresponds to holons, the other corresponds to spinons. This is the first hint of 
the phenomenon of spin charge separation. The full hamiltonian is H = H c + H s 
where, 

H - \^f, 4 tfeT 4 feT -Vf, /ltfc+9t 4 fc +<?T 
n c ~ 2-^, k k-<i~\ k-q^ / j efc A kt ^feT 

kq kq 
kq kq 



k,k ,q 
k,k ,q 



1 especially if np{p f) = np(p I) in other words if the number of up-spins is equal to the 
number of down spins M = M 
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+w a E ««t +4r 9T )(^> g+Gi + 4V 9+Gi ) (is) 

fe,fe ,g 

kq kq 

+ fe fe-«T fe-«T k feT *=T 

/eg /eg 

-£E(4it+^r T )(4 T - 9i +4f- 9i ) 

_ 17 f/l fc 'j- 4- 4tfc+9Tv4 fe T , Atk'-q-Gl-, 

N a 2^ y A k+q] + A k 'l >y A k'-q-Gl +A kt ) 

k,k ,q 

-w a £ i4u + A l>t q ^ Ak J- q+Gl + 4~ q+Gi ) as) 

^1^°" destroys an electron with some momentum and creates another with a dif- 
ferent momentum without changing the spin. This means that this corresponds 
to charge fluctuations. A\° destroys an electron with some momentum and 
creates another with a different momentum, and in addition changes the spin 
projection by ±1. This means that this corresponds to spin fluctuations. If we 

want the luxury of intepreting the objects A\£ as being exact bosons obeying 
the canonical commutation rules : 

[4/,4t V ] = n F (k)(l-n F (k)) 9(n- \k\)9(n \k'\) (20) 

(and all other commutation rules involving any two of these objects are zero) 
then we must pay a price in that the momentum distribution is somewhat com- 
plicated, namely it is appropriately resummcd in order to tame the infrared 
divergences that are known to be ubiquitous in one dimension. Thus the mo- 
mentum distribution is given by, 

< n ka >= npifyF^ka) + (I - n F (k))F 2 (ka) (21) 
Fl (ka) = X - (l + e- 2S «( feff )) (22) 

F 2 {ka) = ^(\-e- 2S< ^) (23) 
S%{ka) = ]T < G \Al k +^A k +^\G > (24) 

S° A (ka) = J2< GolAt^A^JGo > (25) 

Here |G > is the ground state of H = H c + H s . In what follows, we assume 
that M = M , that is, the total number of upspins is equal to the total number 
of downspins. 
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3 Solution Away From Half-Filling 

We focus on the situation away from half-filling where the diagonalisation is sim- 
plest since we may ignore umklapp processes. The holon part of the hamiltonian 
is given by, 

-"c — 2^ fc fe- 9 T fe-?T fe fcT feT 

kq kq 

+E e * 4 fe 44i gi -E^ < +qiAk H ql 

kq kq 

+w a E +4r ,T )(^i 9i +4'f 9i ) (26) 

and the spinon part is given by, 

+ 4 tfeJ - 4 feJ - -Vf, A^ k+ql A k+ql 

+ fe fe-gT fc-«t fe feT *=T 

fcg /eg 

E (<V + 4ft' T )(4?- 9i + 4~ qi ) ( 27 ) 

We use the equation of motion method to solve for the boson propagators. 
Since the details are tedious and the method has already been highlighted in 
our earlier work[l], we shall merely write down the final answers for the boson 
occupation. 

A \ka ( ,+ , A ka l^\_ U2 1 \ ^ P fa iz n) \ A k-q a , A k?qo\ 
A-qAt )A k -qAt)) ^,—^ tc ^ lZn) { - lZn - €k+€k _ q? ^ 

A^<? (t +\ A ka (f] \ U i 1 [ A k-q3i A l-q*\ (2q] 



r \ka &\k(7 i \ \k-qc Ak-qo^ 

p{ %Zn) = ^ l^gl^gj + ^ A k« ,A ] 

{-IZn - e k + £k-q) (lZ n - dk-q + Cfe) 



(30) 



u 2 

e c (q, iz n ) = 1 - -TH p2 (l^ iz «) ( 31 ) 

■'"a 

e s (q, iz n ) = 1 + ^-P(q, iz n ) (32) 

At absolute zero, the discrete sums become integrals and we have to evaluate 
the integral by transforming it into a contour integral in the complex plane such 
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that only zeros of the dielectric function ( and the pole of P ) are included. This 
means that we have to close the contour such that the zeros of the dielectric 
function are found at iz n > 0. Thus we may simplify the expressions as follows. 



qa J 



q ^ w, Nle c { q ,^)£l^P-^q^){-^(q)-^ + ^- q ) 2 

U 2 1 [ A k-qa> A k k -q<j] ,„„s 

Kp-i(q,u; c )£l= Uc e c (<l,")(-"c(q)-e k + e k - g )i ( > 

Here w* > is the zero of P _1 and lo c > is the zero of e c . 

P- I (g,^) = (34) 

e c (q,co c ) = (35) 
Similarly we may evaluate the other correlation function as, 

4V* + )^(*)) = --£ ,i 1 , u j^rrr ^ (36) 

and w s > is the zero of e s . 

e s (q,Lo s ) = Q (37) 

Since we are interested only in the asymptotics we focus on the small \q\ regime. 
In this regime, 

P(g,*)«£ 7 /T^/U (38) 
2tt (uj 2 - <isin 2 (k F )q 2 ) 

Since e c {q,w*) = oo we may ignore this contribution. Let us now enumerate all 
the zeros of the dielectric function. 

/ \ U 2 (4 sin(k F ) q 2 ) 2 

tc(q,u) = 1 - 1 F q > 2 =0 (39) 

(^tt) - 4,sin 2 (k F )q 2 ) 

U (4 sin(k F ) q 2 ) „ , , . 

There is only one zero of e s namely, 

w s (g) = v F |g| (l ) =v 8 \q\ (41) 

where vf — 2 sin{kp)- Thus we may write the spinon hamiltonian in terms of 
the elementary excitations of the interacting system as, 

H s = J2^(q)dl(q)d s (q) (42) 
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and for the holons, 



u c ,i{q) = v F \q\ ( 1 - — ) =v Ci i\q\ 

7TVp 



^0,2(9) = v F \q\ 1 



v c ,2\q\ 



(43) 



(44) 



H c = ^2w Cil (q)dl A {q)d Cll {q) + ^u c , 2 {q)d\ 2 {q)d ca {q) (45) 

It appears that these formulas are valid only for weak coupling. For strong 
coupling we shall have to find some other way. For now we focus only on weak 
coupling. 



U 

A l k -qA t+ ) A k-qa( t )) = 1^2 E 



1 



\Aka A-fka 
l^-k-qcn A k-, 



N 2 



« W P_1 (9.Wcj)jejI"="c,jCc(«,w) (-w c ,j(g) - £fe + £fc- 9 ) 2 



E 



c/ 2 



[/lfe<7 /tt fcCT 1 



4?riV a u C)J - (-VcjM sgn(k)) 2 

U \q\ [A k k % s ,A\%] 



N a -(2it/U)(v s /vf) {-v s \q\ -v F q sgn{k)) 
Therefore for k close to +kp we may write ( k = k F + x ), 

\\ka A k 
j ^ ^k-qa^k-qa ^ T ^ ^ ^ 

9 9 



(46) 
2 (47) 



>+E<4-V^> 



A dq y, U 2 v F 



1 



7^ 2 87f2 W cj («c,j 



+ 



r A dq 

J\x\ q 



u 2 v F 



4n 2 v s (v s + v F ) 2 2 



lLog 



The anomalous exponent is then given by, 

U 2 v F 1 U 2 v F 1 

'-f 2 4tt 2 v c j (v c j +v F ) 2 2tt 2 v s (v s +v F ) 2 



7=E 



' A 

. I x I 
(48) 

(49) 



The velocities of spinons and holons are consistent with the ones obtained by 
traditional bosonization methods. The anomalous exponent is proportional to 
U for small U both in the work by Shulz[5] and our work. The momentum 
distribution is then given by ( ||fc| — k F \ < A), 



n k = i - ^sgn(\k\ - kp) 



\\k\-kp\ 



(50) 
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For a detailed comparison with the results of Shulz[5] we note that Shulz points 
out that the charge stiffness may be written as a power series in U for small 
U as K p = 1 — U/(ttvf) + .... This means that we may read off a formula 
for the anomalous exponent 7 = (K p + l/K p — 2)/4 w U 2 /(2ttvf) 2 . From 
Eq.( 49) we find that since v s ~ v c ,i s=a v c .2 ~ vp, we may conclude that 
7 rs U 2 1 (4:Tr 2 Vp). This agrees exactly with the Shulz result above. In the large 
U- limit, the velocities of the spinons becomes imaginary and we are unable to 
make progress. This means that the cubic and quartic terms that the above 
approach ignores are now relevant. That is, their inclusion alters the value of 
the exponents. At half-filling again there seem to be some difficulties which we 
have been unable to overcome. Perhaps future publications will address these 
issues. 



Momentum Distribution in Two Dimensions 



The main advantage of the sea-boson approach is the ease with which one may 
generalise the above results to more than one dimension by simply promoting 
the wavenumbers to wave- vectors. 



P(q,iz n ) = 



U 2 



N 2 -0 
U i 



E 



P(q,iz n ) 



f Ak<r 4tk<r 1 



e c (q, iz n ) {-iz n - e k + e k - q ) 2 



E 



[-4 



ka 



A 



fkff 



k— q<T> k— q<x 



N a -i/3 ^ e s (q, iz n ) (-iz n - e k + e k - q ) 2 



[A 



ka 4^ ko " 
k— q<r' k— qcr 



(-iz„ 



(k + Ck-q) 

U 2 



r tk-qa ^fk-qcri 



(iZn - £k-q + £k) 



e c (q,zz„) = 1 - ^P 2 (q,iz n ) 
e s (q,zz„) = 1 + —P(q,iz n ) 



(51) 

(52) 

(53) 

(54) 
(55) 



For evaluating the low energy, long-wavelength limit of the above expressions, 
we may ignore the non-zero lattice spacing and instead focus on the continuum 
limit (which is valid for Uf <C 1 since we have set a = 1) where the expressions 
for the RPA-polarization have already been derived by Stern[12](here 2m = 1). 



Re[P](q,w) 



mk F A |q| 
" 27r|q| ^k F 



(K 

\2k F 



mui 
k F \q\ 



-c+ 



( Jqj_ mu> 
\2k F k F \q\ 
(56) 
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TTIUJ 



2kF &F|q| 



where, 



C± = sgn 



|q| mw 



C± = 0, D± = l 



Jte[P](q,a;)«-^{l 



D± = 0, 

2kF &F|q| 

u> 



|q| _ to 



fejr|q| 
< 1 



2?r y/w 2 - wf,q 2 

UJ 



} 



= 1 + ^{1-- 

2?r \/u 2 - vpq 2 



} 



rjoL _ 










(57) 


> 1 


(58) 




(59) 




(60) 




(61) 



for uj > vf\<i\. 



Veff = Vp 



m 



d_ 

div 



Amv 2 F q 2 
i?e[P](q,^- 2 ^ 2 _ 4q2)3/2 



(62) 



(63) 



This means, 



u 2 



\Ako /itkcr , 



N 2 P- 1 (q,Wc)£L=a; c Cc(q,w) (-Wc(q) - £k + e k - q ) 2 



4 k -V(* + )^k^w) = 



(64) 
(65) 



and w c (q) = i> e //|q|. For a triangular lattice the energy dispersion may be 
written as, = —2 [cos(k x ) + 2 cos{k x /2) cos(V3k y /2)]. Let us assume that 
kFd <C 1 (since we have set a = 1, we must have kp <C 1) in which case the 
Fermi surface is a circle, 



1 



T/lkcr jfkCT i 

l^k-qo-' ^k-qcrJ 



k- qCT l J k- q „U/ 2A|_ ci?e [P] (q , w)( -c c (q)- ek + ek _ q ) 2 
M " (-«e//|q| -k 2 + (k-q) 2 ) 2 

T|q|(«e// - 4) 3/2 0(A* - |k - q|) 0(|k| - k F ) 



Amvp 



(-« e// |q|-2k.q + q 2 ) 2 
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*M£ff ~ v lf' 2 0(k 2 F - k 2 + 2|k||q| C0S (g) - q 2 ) 0(|k| - k F ) 
Amvp (-u e //|q| - 2|k||q|cos(a) + q 2 ) 2 



(66) 



pTT/2 r-OC 

E 4-V(* + KVW =2 da dq 
„ x 'Jo Jo 



t/2 i-oo (^2 ,,2)3/2 g (fc 2 _ fc2 + 2fcFgcog(a) _ g 2) g ( | k | _ fc^ 



27TW 2 , 



K 2 /f-4) 3/2 



- 2- 



ArcTan 


V e ff-VF 


.V v °ff- V F- 





\J v hf- v 



Zp = exp 



TTVi 



4 ArcTan 



V ef f~ V F 
\f V lft~ V F 



"eff 



v hf - v f 



For U — > +oo, we have, 



In this limit, 



Zp(U — > +oo) w exp 









2w F 










8^ 





For [/ — > +0, we have, 



w e // &v F [1 



?7 2 \ 
32tt 2 J 



Z F (C/ -» 0+) S3 1 - U 3 



64tt 4 



(-w e // - 2fc F cos(a) + g) 2 
(67) 

(68) 

(69) 
(70) 
(71) 

(72) 
(73) 



The above results are valid when it is legitimate to approximate the shape 
of the Fermi surface by a circle. Strictly speaking therefore it is valid only 
for fcp <C 1. Also the U — ► +oo result must be taken with a grain of salt 
since the corresponding Id case has not worked out. Just as the anomalous 
exponent saturates to a value of 1/8 according to Bethe ansatz and also hopefully 
according to the sea-boson method provided we include the anharmonic terms, 
we may suspect that the quasiparticle residue also saturates (to a value close to 
unity) when these anharmonic terms are included. This is completely consistent 
with the results of Castro-Neto and Fradkin (who predicted a value of 0.78 I 
think for the smallest value of the quasiparticle residue). However it is not clear 
to the author if their results were due to a formalism more powerful than the 
sea-boson method, or because of a clever use of intuition. In particular one 
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has to see if the traditional bosonization method can predict the exponent of 
1/8 in Id, without relying on the Bcthe ansatz results. The small U results 
are probably reliable though one has to check this several times more. The 
quasiparticle residue seems to deviate from unity by a term proportional to U 3 
rather than U 2 . 



5 Full Propagator 

To compute say the density of states, we have to compute the dynamical one- 
particle Green function. For this we have to complete the construction of the 
field operator started in the Appendix. In particular we have to find a formula 
for R(pa). 

c pCT = e^-i H^ +< " 1+A p+<-i) T " (p) R( pa ) (74) 

c pCT ee e^-'i H^+^Xj Tq(p) R*( pa ) (75) 

The main problem with these formulas is that it is not possible to use the 
Bakcr-Hausdorff theorem since the commutator of the two exponents does not 
commute with either even at equal times. Therefore it is not possible to sys- 
tematically deduce R by making contact with the free theory. However we may 
suspect that perhaps, 

R{pa) =n F {p) e tN ° « (pff) (76) 

as usual tif(p) = 0(kp — |p|). Here £(pcr) is an arbitrary non-constant function 
(c- number) of the arguments such that, 

e -iN° {WgiW «pV) M Sp p , 6„j (77) 

Therefore, 

c pv (*)c pff (0)) = /e^A p " ) >e^A A *° + A p+^J T ^ p ' R *( p a ) ^(p^) 

(78) 

For reasons already mentioned this is not possible to simplify further. 



6 Conclusions 

We have computed the asymptotic features of the momentum distribution of 
the one-band Hubbard model in one and two spatial dimensions. In the one 
dimensional case, we obtain a Luttinger liquid with an anomalous exponent 
identical to the one obtained by Shulz[5] from Bcthe ansatz. This is valid for 
weak coupling only. We are unable to extend the results to strong coupling. In 
the two dimensional case, we obtain a Landau Fermi liquid with quasiparticle 
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residue deviating from unity by a term proportional to the cube of U. Using 
the remarks of Shulz, we may conclude that the two dimensional case is always 
a Landau Fermi liquid. 

7 Appendix : The Field Operator 

In the past, we have focussed on the density phase variable ansatz (DPVA) as 
the main way of ascertaining the connection between the field operator and the 
sea-displacements. This is done in a two step process. The first step involves 
expressing the fields in terms of currents and densities. Next we relate the 
currents and densities to the sea-displacements. Unfortunately this program 
has faced difficulties since it appears now that the phase functional is extremely 
nontrivial thus making the approach practically useless. Here we instead focus 
on trying to pin down the field in terms of the sea-displacements directly. This 
approach is better since we now have a definition for the sea-displacements 
in terms of the Fermi fields. Consider the commutators(we have dropped the 
square roots of the number operators in the denominator), 

[c pCT , A kCTl (qo- 2 )] w Vr- q /2, k nF(p)(l ~~ n HP + q)) T q(p) c p^ <Wi ( 79 ) 

[cp^Lnfa ^)] « S pM+cl/2 n F (p - q)(l - n F (p)) T_ q (p)c p <5 CT;IT2 (80) 

where T q (p) = ea;p(q.V p ). We may now exponentiate the commutation rules 
in Eq.( 79) and Eq.( 80) and arrive at, 

Cva = e E q< / H: +q/2 ^')+ v,/-' t ^ R{pa) (si) 

c pCT = e E ^' ^ R * {p(T) (82 ) 

Here R(pa) is a c-number function that commutes with everything. To verify 
the commutation rules we proceed as follows. The exponential is interpreted as 
a power series expansion. In the expansion, the translation operator translates 
all the p's to its right by an amount q and not just the p in R(pa). 

[c pff ,A kffl ( q( T 2 )] = e ^A- A U^M°')+\ +Ci/2 A-<v)) T ^ p ^ kCTl (q<7 2 ) R(pa) 

_ Akffi(q(T2)e E q .'(-^w-^')+' 4 P +q / 2 .'(-^)) T ^ R {pa) 

= -X^[^P+q/2 ( r( C l Cr ')' j4 k ffl (qcr2)]Tq(p) C prJ 
a 

= ^k,p+ q /2^ llCT n F (k - q/2)(l - n F (k + q/2)) T q (p) c pa (83) 
[c P .,4 CTl (q^)] =e^A- A l + ,^')+\ + ^(-^) T " (p) 4 CTi (q CT2 ) R(pa) 
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-4 CT1 (q^)e^-'(-<+o/- (q<7 ' )+A p+<./-' ( - q<7) ) ^ R{pa ) 
= 5Z[ A P- q /2-' fa* 7 )' (q^)]T_ q ( P ) Cpff 

/ 

= <5k,p-q/2^,a 4 n F (k - q/2)(l - n F (k + q/2)) T_ q (p) c pCT (84) 

In the sea-boson method, unlike in the usual bosonization in Id, the momentum 
space description is easier than the real space description. Now we have to 
compute R(pcr) by making contact with the free theory. This is proving to 
be harder than expected since the commutator of the exponent with the other 
exponent does not commute with either exponent, hence we may not use the 
truncated version of the Baker-Hausdorff theorem. If one insists on computing 
the full dynamical propagator which is needed for obtaining the density of states 
as a function of the energy, then perhaps a simple route may suffice. It goes by 
the name 'serendipitous surmise'. 



7.1 A Serendipitous Surmise 

Many years ago the author had a conversation with then the student now Prof. 
A.H. Castro-Neto where the latter suggested that maybe the field operator is 
simply given by, 

V>(x)«e" < ^ e,qXXq v /^ (85) 

where by definition Xq = iq • j(— q)/(q 2 iV°) and p° = N°/V. Later he realised 
that if we choose LY q , X /] = as is in fact mandatory 2 , then fermion commu- 
tation rules are not obeyed. However one may take solace in the fact that at 
least one commutation rule does come out right namely [^(x), p q ] = e lq x -0(x). 
A refinement over this ansatz was attempted in our earlier work[l] by introduc- 
ing an additional phase functional of the density linear in the density and this 
was also inadequate since by now the author knows that <& there when com- 
puted was imaginary when it was postulated to be real. A compromise was also 
suggested that involved multiplying and dividing by the free propagator and 
using the exact version in the numerator and the bosonized free propagator in 
the denominator. This trick though repugnant to most, gives us an anomalous 
exponent of the Luttinger liquid as we shall see below. However, this anomalous 
exponent is off by a factor of two from the exact one obtained by Mattis and 
Lieb. From Eq.( 85) we may write, 

G(x-x') = (^( x ')^( x )\ r, / e ^o( eq ~ E<1 ) X A p ° 



■£ q (l- COS [q.(x-x')])(X q X_ q } 



P v (86) 



2 Since strictly speaking it is the conjugate to p as defined by the line integral of the ratio 
of the current and density that enters; these commute amongst themselves. 
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Again using the trick outlined in our earlier work[l] we may write, 

G(x-x) = Go(x-x) e -£,(i-^qW)])[<*,*-,M*,*-,> ] (87) 

Here G (x — x ) is the propagator obtained from elementary considerations. In 
one dimension, we may see that (Y q Y_ q ) « kpS(q)/(q 2 N°). The structure 
factor So(q) = \q\/(2kp) for the interacting case. For the interacting case we 
have, S(q) = {v F /v eff ) S (q). 

f°° J„ l-ch.^-.')] I v f A 

G(x - x) = Go(x - x) e 9 ra V v °" V 

where 7 = — 5 • This exponent is exactly one half of the exponent obtained 
by Mattis and Lieb. What is even worse is, we have shown in an earlier preprint 
that when applied to the X-ray edge problem, we obtain the well-known results 
of Mahan in one dimension but not in higher dimensions. Thus it would appear 
that there is something amiss in the expression for the field operator. Neverthe- 
less it may be a quick and easy way of getting the density of states that gives 
the right qualitative physics if not the right exponents. 
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